Nonequilibrium phase transition due to communities isolation 
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We introduce a simple model of a growing system with m competing communities. The model 
corresponds to the phenomenon of defeats suffered by social groups living in isolation. A nonequi- 
librium phase transition is observed when at critical time tc the first isolated cluster occurs. In the 
one-dimensional system the volume of the new phase, i.e. the number of the isolated individuals, 
increases with time as Z r^ t^. For a large number of possible communities the critical density of 
filled space equals to pc = {m/NY'^ where A'^ is the system size. A similar transition is observed 
for Erdos-Renyi random graphs and Barabasi-Albert scale-free networks. Analytic results are in 
agreement with numerical simulations. 

PACS numbers: 05.50.-|-q, 89.75.Hc, 02.50.-r 



Recently physicists working on modeling of social phe- 
nomena are frequently touching the idea of dissemina- 
tion and competition - especially in the case of language 
[]|, culture y and opinions [3|. The key subject can be 
posed as follows: how does the node internal variable 
change when it is influenced by others? The issue has 
occurred especially in Axelrod model of culture dissemi- 
nation VM or many sociophysics systems such as: Sznajd 
model y, voter model Q, majority rule voting @j_social 
impact model '?| or bounded confidence models [8|. At 
this moment one should stress that the above mentioned 
models are bound to explore social effects of diffusion or 
adoption of node states. Such processes are usually run- 
ning in accordance with the following scheme: one takes 
the state of its neighbor, provided that a set of rules is 
fulfilled. 

There is however also another, qualitatively different 
phenomenon - isolation of surrounded social groups and 
resulted extinction of their members due to lack of com- 
munication with other groups of the same specie. In fact 
this issue should be regarded as equally important as the 
previously mentioned dissemination or migration effects. 
The isolation and consequently the lack of communica- 
tion among the groups belonging to the same community 
(and vice versa - no communication causing the isolation) 
might lead to severe disturbances in the society. One of 
them can be racial isolation (segregation) that can cause 
serious social problems [9|. In fact the phenomenon of 
residential segregation has been studied in several phys- 
ical papers |lOl| that in part follow the famous work of 
Schelling |ll| or other Ising-like approaches. In other 
situations the lack of social contact effects in increased 
mortality of seriously ill patients as compared to those 



that are not isolated [l4|. Finally, recent research [1^ 



shows that Americans suffer from social isolation due to 
dramatic decrease of number of discussion partners even 
with those that they share the closest relationship. Given 
the fact that such social phenomena as elections [IJ] or 
war [la ] are currently being examined using methods of 



statistical physics a quantitative model of social isola- 
tion could be a useful tool to predict blocking of voting 
districts or trapping of hostile troops during a wartime. 
The key idea of our work can be presented in form 
of two questions: (i) what happens if instead of species 
spreading the interaction effects in species isolation and 
extinction! (ii) what are consequences of the fact that 
sometimes a small group is capable to surround and de- 
feat a larger one? The first point comes as an effect of 
the observation that a group of people that is suddenly 
surrounded by people from opposite groups is often de- 
feated by enemies or opponents. It occurs because the 
surrounded group is isolated and is not able to communi- 
cate with other group members who could support them, 
e.g. provide a military backup (it was a common case 
during many wars). To justify the second question one 
can think of Chinese game Go where, in some condition, 
one player can block the opponent using the amount of 
stones which is less than those that are just being sur- 
rounded. 

t = 2 






FIG. 1: (Color online) An example of evolution in the chain 
consisting of 8 nodes. Open circles are empty sites, black 
and gray circles correspond to different communities. Isolated 
nodes are marked with a cross. 



In this paper we impose these dynamical rules onto 
various regular and random networks. We start with a 
simple chain where sites are being filled with individu- 
als belonging to two different species. Then we extend 
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FIG. 2: (Color online) (a) Number of isolated nodes (Z, dotted and dashed lines) and not isolated nodes of each specie 
{{t — Z)/2, solid lines) versus time for different chain sizes A'' (gray dotted - fO^, black dotted - 10'', gray dashed - 10^ and black 
dashed 10®). (b) Number of isolated nodes (Z, filled symbols) and not isolated nodes of each specie {{t — Z)/m, open symbols) 
versus time for different number of species m (circles - m = 4, squares - m, = 16 and triangles - m — 64). All simulations are 
for A'^ = 10^ and the lines come form the solution of Eq. ([5} 



the model to a case of m different species and finally 
we consider the Erdos-Renyi graphs and Barabasi- Albert 
scale-free networks. 

Let us consider a chain of N initially unoccupied nodes. 
In each time step one empty node is chosen randomly. 
Then, an internal variable (t) or (|) for this node is ran- 
domly selected. Both possibilities correspond to differ- 
ent species or communities and are drawn with the same 
probability. If a cluster of n identical filled nodes (e.g. 
tlTT) is surrounded by individuals belonging to other 
community (e.g. iTTTTi)i the nodes in the surrounded 
cluster are treated as extincted and can no longer inter- 
act with the rest of the chain, i.e. they will be not able to 
surround other clusters. The procedure is held until the 
chain is full, which happens in the iVth time step. An 
example of complete evolution of the system is presented 
in Fig. m 

Our main points of interest are: (i) the critical time tc 
when the first isolated cluster appears (ii) the number of 
isolated nodes for t > tc- 

Figure [2^ shows the number of isolated nodes Z and 
the number of not isolated nodes of both species {t — Z)/2 
as function of time for four different chain sizes: 10'^, 10'*, 
10^ and 10^. In each case the number of isolated sites 
follows a power law Z ~ t" with a exponent close to 3.0 
{a = 3.09 for N = 10^, a = 3.06 for N = 10'', a = 3.05 
for N = 10^ and a = 3.04 for N = 10*^). 

The plots indicate that in this system we observe a 
nonequilibrium phase transition - after reaching a certain 
time of the evolution (after filling a specific number of 
nodes) a new phase emerges due to the occurrence of 
the first isolated cluster. The volume of this phase can 
be treated as the system order parameter. It grows up 
when one runs above the critical time tc. Moreover it 



can be seen in Fig. [3^ that data for different chain sizes 
collapses onto one curve after rescaling both Z and t axis 
by system size N. Figure [3)d shows that the critical time 
tc of the first isolated node appearance grows with system 
size as tc ~ A^" where f3 = 0.664 ± 0.001. 

In order to obtain the average number of isolated 
nodes, we have to sum all different possibilities of a clus- 
ter to become isolated. A single isolated site emerges ei- 
ther as an effect of a combination f jt or |t| in which the 
middle node is turned into isolated one. To express the 
total number of such nodes in the system (Zi) we need to 
multiply the probability of the sum of those combinations 
by the number of such possibilities, that is A^ — 2. Sim- 
ilarly the number of isolated sites coming from isolated 
clusters of size n is: 



Zn=n{N-n-l) Pr2(t) J] P<i) + Pr'(i) H P''(t) 

1=1 1=1 

(1) 
where n = 1, 2, . . .. As the examined system is symmetric 
(i.e. Pr(t) ~ Pr(|)), taking into account that at time t 
there are already Z isolated nodes the average probability 
of finding a certain specie at time t is {t — Z)/{2N). Since 
Z = X)i=i ^i ^"2 obtain after short algebra: 



Z = 2 ^ (n-2)(A^-n+l) 



n=3 
^3 , 



t- z 

2N 



(2) 



or Z ^ {t - Zf/{2N - t + Zf. Solving this equation 
leads to 



1 



Zr^- \Ur - 4 + (8 + IQtr - tl)Ur 

6 '■ 



(3) 



with Ur = (v/3(16 - 24t^ + mtl - 2tl) - 80 -K 84i,. - 
24t2)i/3 where t^ = t/N and Z,. = Z/N. The formula 
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FIG. 3: (Color online) (a) Data collapse for rescaled number of isolated nodes (Z/N) versus the rescaled time (t/N) observed 
for three different data sets - N = 10^ (circles), A^ — 10^ (triangles) and A*' — 10® (squares). The curve (hardly visible) is 
obtained from Eq. @. (b) The critical value of time t^ for m = 2 versus the chain size A'^. The line is Eq. Q while triangles 
are numerical simulations, (c) The critical value of time tc versus the number of species m for different chain sizes A'^ — 10'^ 
(circles), A' — 10^ (triangles) and A'^ — 10^ (squares). Lines come from Eq. ((SJ. 



is universal for any value of chain size - all data should 
collapse on this curve, as it can be seen in Fig. [3^. If 
2' <C t <C iV, what is acceptable for the most part of 
the evolution, then Eq. ^ leads to Z « t^/(47V2) i.e. 
the number of isolated nodes should increase as t^ . This 
fact is in agreement vi^ith the numerical experiment. This 
approximated formula can be also used to calculate the 
critical time tc at which the first isolated node appears 
in the chain. Putting Z = 1 we get a simple expression 
for the critical time 



tc - i2Nf' . 



(4) 



This result is consistent with the value of the /3 parameter 
observed in the numerical data. 

We can easily extend the previously described model 
of two competing species onto a case where the number 
of species is tti > 2. Similarly to the two-species case, in 
each time step a type of specie is drawn from the uniform 
distribution (l,?n) and placed in a random, unoccupied 
place in the chain. The isolated nodes are formed from a 
cluster of identical species surrounded by other identical 
species. 

An example of the evolution of the extended model is 
presented in Fig. [^h- Like in the case of m = 2, the 
number of isolated nodes follows a power-law Z ^ V 
with 7 exponent close to 3 (7 = 3.05 for m = 4, 7 = 3.00 
for m — \% and 7 = 3.01 for m — 64). 

The analytical approach in the case rn > 2 is identical 
to the case to = 2 except for two things. First, there 
are m different species which can be isolated. Each of 
those m species can be isolated in ttj — 1 ways, therefore 
instead of factor 2 we should put m(ra — 1) in front of 
each equation in the set of equations H]). Second, larger 
number of species results in the change of probability of 
finding a specific specie - in the extended model it is equal 
to (i - Z)/{mN). Thus Eq. ^ has now the form 



Zn = n ■ m{m — 1){N — n — 1) 



niN 



n+2 



where n = 1,2,... . Following an identical algebra as 
in the case of two-species model, we arrive at a self- 
consistent equation for the number of isolated nodes 
Z = (m - l)(i - Zf/{mN - t + Zf which is exactly 
algebraically solvable. The solution fits to the numerical 
data quite well (see Fig. ^p) and as before one can ap- 
proximate it with Z « {m — l)t'^ / (jnN)'^ . This proves 
that the increase of isolated nodes follows the same rule 
as in the two-specie case, i.e. t^. The critical time t^ at 
which the first isolated node appears is 

tc^(^) N^l\ (6) 

\TO — 1/ 

which, once again, is consistent with the numerical data 
(see Fig. [St). 

The form of Eq. ([6]) gives us the opportunity to spot 
the interplay between the only two parameters of the 
model - the length of the chain N and the number of 
species m. If to ^ 1 the Eq. ^ can be rewritten in a 
form of tc = {mN'^y/^, what leads to the critical density 
of filled nodes pc = {m/Ny^^. The obvious conclusion 
from this relation is that when the chain becomes larger 
the critical density gets smaller and in the thermody- 
namical limit vanishes completely. Then, if we would 
like to maintain a constant value of pc we should require 
m/N — const. In other words, it is possible to prevent 
the convergence of critical density to zero by making the 
number of species proportional to the chain's length. 

The critical time t^ for to, = 2 can be also found in the 
case of two- and three-dimensional cubic lattices and, 
what is more important, for any random network char- 
acterized with a specific degree probability distribution 
p{k) (fc is the number of links of a given node). The 
general fornuila for number of single isolated nodes is 



Zi = 27V ^ p(fc)x'=+\ 



(7) 



fe=0 



(5) 



where x = t/{2N). In order to obtain the critical time, 
we require that Zi — I and solve this equation for t. In 



case of regular lattices we have tc = {2N)^/^'^~^^\ where 
z is number of neighbors. We have calculated the crit- 
ical time for two most popular types of complex net- 
works: Erdos-Renyi random graphs [l6| and Barabasi- 
Albert scale-free networks |17|. In case of ER graphs 
characterized by degree distribution p{k) = e^^*^^ (k) /k\ 
the critical time Ieb. can be expressed as 
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27V, 
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(k) 



2N 



(8) 



where W{x) is Lambert W-function. For BA network 
(degree distribution p{k) = \{k) k~^) we get 



tsA 



'4)^\2NY 



(fc) 



(9) 



The obvious condition for avoiding an isolated node in the 
system is tc > N. It leads to the following inequalities 
for different networks: N > 2^ for regular lattices, N > 



(fc>- 



" (fc) for BA networks and N > e^*^^/^ for ER graphs. 
The above described results are shown in Fig. [H 




FIG. 4: (Color online) Logarithmic plot of the critical time tc 
for Erdos-Renyi graphs (empty symbols) and Barabasi- Albert 
networks (full symbols) for different networks sizes: A'^ — 10"' 
(circles), A'^ = 10'* (triangles) and A'' — 10^ (squares). Sym- 
bols are numerical simulations and lines come from Eqs. (|8} 
and (|9|. The inset shows tc in linear scale for A'^ — 10^ in 
case of ER graph (triangles) and BA network (squares) - tc is 
rescaled by a factor of 10000. 

Conclusions - In this work we proposed a simple ap- 
proach to model communities isolation in growing soci- 
eties. The numerical simulations, fully supported by an- 
alytical approach show that a critical time tc a nonequi- 
librium phase transition takes place and a new phase 
consisting of surrounded clusters emerges. In the case 
of one-dimensional system the number of isolated nodes 
rises with time as a power-law with exponent 7 = 3. The 
scaling is universal i.e. it depends neither on the chain's 



length N nor on the number m of possible species. An 
analytic form for the critical time tc is found and for 
large m this time scales as tc — {mN^y/'^. The phe- 
nomenon has been also observed for higher dimensional 
systems as well as for Erdos-Renyi random graphs and 
Barabasi- Albert scale-free networks. 
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